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Abstract 

We consider the capability of the e + e~ linear collider (which is recently called as the 
International Linear Collider, or ILC) for studying the properties of the heavy Higgs 
bosons in the supersymmetric standard model. We pay special attention to the large 
tan (3 region which is motivated, in particular, by explaining the dark-matter density 
of the universe (i.e., so-called "rapid-annihilation funnels"). We perform a systematic 
analysis to estimate expected uncertainties in the masses and widths of the heavy Higgs 
bosons assuming an energy and integrated luminosity of y/s = 1 TeV and L = 1 ab . 
We also discuss its implication to the reconstruction of the dark-matter density of the 
universe. 



1 Introduction 



An e + e~ linear collider, the International Linear Collier (ILC), has been discussed as a new 
experimental tool to study the physics at the electroweak scale and beyond [H El EH E| • The 
primary purpose of the ILC is to clarify the physics of the electroweak symmetry breaking 
and physics beyond the standard model; in order to solve the naturalness problem in the 
standard model, it is expected that some new physics shows up at the energy scale of 100 GeV 
— 1 TeV. Discovery and, in particular, precise studies of such new physics should become 
very important once the ILC will be built. Thus, at the current stage, we should study the 
potential of the ILC for such purposes. 

Among various possibilities, supersymmetry (SUSY) is one of the prominent candidates 
of the new physics beyond the standard model. Accordingly, superpartners of the standard- 
model particles are regarded as significant targets of the ILC. As well as the superparticles, 
however, we should also study another class of new particles which show up in the supersym- 
metric models, that is, the heavy Higgs bosons. Since the minimal supersymmetric standard 
model (MSSM) contains two Higgs doublets, i.e., up- and down-type Higgs bosons, CP even 
and odd neutral Higgs bosons (H and A) as well as the charged Higgses exist in the 
physical spectrum as well as the standard-model-like Higgs h. Detailed study of the heavy 
Higgses provides information about some of the parameters in the MSSM like, for example, 
tan j3 El HI El El Hill • Thus, the study of the heavy Higgses will be important for the study 
of the Higgs potential and electroweak symmetry breaking. 

The study of the heavey Higgses is important both from the point of view of understand- 
ing electroweak symmetry breaking as well as for its deep impacts on cosmology. It has long 
been recognized that the dark matter in the universe, whose origin cannot be explained in 
the framework of the minimal standard model, may be explained in the MSSM. With the 
i?-parity conservation, the lightest superparticle (LSP) is stable and hence, if its interaction 
is weak enough, LSP can be a viable candidate for the cold dark matter. Importantly, in 
large fraction of the parameter space, the lightest neutralino becomes the LSP which can be 
the cold dark matter in this case. In the past, the relic density of the lightest neutralino has 
been intensively studied HJJ [T21 H31 H3J [IHl UHl HTJ HHJ [TU] . 

It should be noticed that, in order to realize the LSP dark matter, there is one caveat; in 
many cases, the thermal relic density of the LSP becomes larger than the currently measured 
dark matter density. In particular, the density parameter of the dark matter is now well 
constrained by the WMAP [HI I2T] 

fi (WMA P)/l2 = o.113±8;88|, (1.1) 

and it is not automatic to realize such density parameter with the LSP. It is often the case 
that some mechanism is needed to enhance the annihilation of the LSP to realize the LSP 
dark matter. Detailed study shows that the relic density of the LSP can become consistent 
with the WMAP value with some particular annihilation process. In the so-called mSUGRA 
models, the parameter regions where the density parameter of the LSP becomes consistent 
with the WMAP value are classified as "bulk region," "coannihilation region," "focus-point 
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region," "rapid-annihilation funnels," although the annihilation mechanisms used in those 
regions work in more general framework. 

Once the supersymmetry is discovered, one of the challenge will be to determine the 
thermal relic density of the LSP and to see if the LSP dark matter is plausible |2*2*| EHJ 
EH ESI ESI EH EH EHj- If successful, it will give us deeper understandings of our universe 
up to the temperature of O(10 GeV). For this purpose, it is necessary to measure various 
parameters in the MSSM for the calculation of the thermal relic density of the LSP. Such 
studies have been also performed for the case of the LHC [3U1 Unj ED E21 ESI Ell- Most of 
those studies have, however, assumed very simple model of supersymmetry breaking, like 
the mSUGRA model, in estimating the expected uncertainties in the reconstructed value of 
the thermal relic density of the LSP. We believe that, once the superparticles are found, it 
is necessary to calculate the relic density in a model-independent way. Only with the LHC, 
such a study seems difficult in many cases. 

With the ILC, on the contrary, precise determination of the MSSM parameters will be 
possible without assuming any particular model of SUSY breaking. Thus, in this case, the 
ILC will be able to help reconstructing the relic density of the LSP in many cases. For 
example, in the focus-point case, in it was shown that the ILC can provide useful and 
significant information for the calculation of the thermal relic density of the LSP. 

Here, we pay particular attention to another case, the rapid-annihilation funnels. In this 
case, pair annihilation of the LSP is dominated by the diagram with the s-channel exchange 
of the CP-odd neutral Higgs boson. In particular, in order to enhance the pair annihilation 
of the lightest neutralino, tan j3 parameter becomes large in this case. If the LSP dark matter 
is realized in the rapid-annihilation funnels, detailed study of the heavy Higgs sector in the 
large tan j3 case is necessary to calculate the annihilation cross section of the LSP. 

Strongly motivated by the dark matter physics in the rapid-annihilation funnels (as well 
as the study of the Higgs potential), in this paper, we consider the strategy for studying 
the properties of the heavy Higgs bosons at the ILC assuming an energy and integrated 
luminosity of y/s = 1 TeV and L = 1 ab _1 . We concentrate on the case where tan/3 
parameter is large and see how and how well we can constrain the properties of the heavy 
Higgs bosons. As we will discuss in the following sections, interactions of the heavy Higgses 
are well parameterized by five parameters. Using pair production processes of the neutral 
and charged Higgs bosons, these parameters can be constrained. If information from the 
neutral and charged Higgs production processes is combined, it will help improving the 
determination of properties of the heavy Higgses. In this paper, we particularly consider 
how to combine information from various processes and how well the resultant constraint 
can be, using systematic consideration of the production and decay processes. Then, we also 
study how the results can be used to reconstruct the thermal relic density of the LSP. 

The organization of this paper is as follows. In Section El we briefly review the properties 
of the heavy Higgs bosons. In Section El we discuss how and how well we can study the 
properties of the heavy Higgses at the ILC. In particular, we estimate the expected errors 
in the observed values of masses and widths of the heavy Higgses. Then, in Section we 
discuss the implication of the ILC study of the heavy Higgses to the reconstruction of the 
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dark matter density. Section is devoted to the conclusions and discussion. 



2 Higgs Sector 

We start by discussing the Higgs bosons in the MSSM and by summarizing the relevant 
parameters for our analysis. As we mentioned, the MSSM contains two Higgs bosons H u 
(1,2, ~) and H d (1,2*,—-), where we have shown their quantum numbers of SU(3)c X 
SU (2) i x U(l)y gauge group in the parenthesis. Both of these Higgs bosons acquire vacuum 
expectation values, which are parameterized by 

(H° u )=vsin(3, (H° d )=vcosf3, (2.1) 

where v ~ 174 GeV, and the angle (3 is a free parameter. After the electroweak symmetry 
breaking, it is usually the case that one linear combination of H u and H d becomes the 
standard-model-like Higgs $sm while the other combination $ contains physical heavy Higgs, 
which are given by* 1 




cos (3 sin (3 
- sin (3 cos (3 




(2.2) 



Using the mass eigenstates H, A, and H ± , the heavy Higgs doublet $ = ($°, $ ) is expressed 
as 

$° = -^=(H + iA), <5> ± = H ± . (2.3) 
v2 

In general two Higgs models, masses of H, A, and H are unconstrained. In the MSSM, 
however, their masses, tuh, tua, and ra#±, are quite degenerate within a few GeV or so. 
For the neutral Higgs sector, this fact makes the study difficult. In our following study, we 
do not adopt the theoretical relation among these masses but we consider how well we can 
determine these three masses (in particular, ran and tua) at the ILC. 

The productions of the heavy Higgses at the ILC are dominated by the gauge-boson 
exchange diagrams. In Fig. ^ we show the cross section for the processes e + e~ — ► AH and 
e + e~ — > H + H~ for several cases. (We will give the formula for the cross section in the 
following section; see Eqs. ()3.5|) and ()3.19|h ) Here and hereafter, in our numerical study, we 
assume un-polarized electron beam. 

Decay processes of the heavy Higgses have more model-dependence. In the rapid- 
annihilation funnels, heavy Higgses mostly decay into the third generation quarks and lep- 
tons. We write the relevant interaction as* 2 

Yukawa = Vb^Qs^R + Vt®* icr 2 Q 3 t R + y T $L 3 T R + li.c., (2.4) 



^Strictly speaking, mixing in the CP-even neutral Higgs sector differs from those of other Higgses, and 
conventionally the mixing in the CP-even neutral Higgses are parameterized by the parameter a |35|. When 
the heavy Higgses are much heavier than the lightest Higgs (i.e., in the so-called "decoupling limit"), however, 
it is expected that sin(/3 — a) — + 1 and Eq. Q2.2|l holds with good accuracy. In this paper, we only consider 
such a case. 

# 2 Notice that the Yukawa coupling constants yb, yt, and y T are different from those in the superpotential. 
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Figure 1: Total cross sections for the heavy Higgs production processes as functions of the 
center of mass energy Eqm = \fs. The solid lines are for the process e + e~ — ► AH with 
rriA = viiu = 300 GeV, 350 GeV, 400 GeV, and 450 GeV from above, while the dotted lines 
are for e+e~ -> H+H- with m H ± = 300 GeV, 350 GeV, 400 GeV, and 450 GeV from above. 



where Q 3 (3,2,|), b R (3*,1,§), t c R (3*,1,-|), L 3 (l,2,-±), and r c R (1,1,1) are third- 
generation quarks and leptons. In the parameter region we are interested in, the (funda- 
mental) Yukawa coupling constants in the superpotential are fairly large. So, we expect that 
the radiative corrections do not significantly affect the relations between those with Yukawa 
coupling constants of the heavy Higgses. If we try to extract the information about these 
Yukawa coupling constants from the observed fermion masses, however, radiative correction 
may become important. At the tree level, the Yukawa coupling constants are related to the 
masses of the bottom quark, top quark and r as 

(tree) _ m fe tan/3 (tree) _ m t cot p (troe) _ m T tan j3 
Vb ~ j Vt — > V T ~ • K^-v) 

V V V 

Thus, once the "bottom-quark mass" is well determined by the study of, for example, the 
lightest Higgs boson h, then the Yukawa coupling constant for r can be also well determined 
at the tree level. It has been pointed out, however, that y& receives sizable radiative correction 
from supersymmetric loop diagrams in particular at the large tan /3 region [SHI EZ] • Thus, 
in determining the interaction of the heavy Higgses, it is dangerous to adopt the tree level 
relations. 

In our analysis, we assume that several qualitative features of the large tan (3 region 
hold; in particular, we consider the case where the heavy Higgses decay into third generation 
quarks and leptons because of large y& and y T . We do not, however, rely on the tree level 
relations among these Yukawa coupling constants but regard all the Yukawa couplings (in 
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particular, y b and y T ) to be free parameters which should be experimentally determined. 

In the large tan (3 region, y t is extremely suppressed and, hereafter, we neglect its effects. 
Then, the neutral Higgses A and H decay mostly into the bb or rf final states.* 3 The decay 
rates of the heavy Higgses are given by 



A^bb 



H^bb 



r, 

^H ± ^bt 



Jk 

167T 
1 

167T 

Ik 

167T 
1 

167T 

Ik 

167T 
1 

167T 



y b m A \1 + l\ bb ) , 



yr m A, 



9t v H (i + Air») 



y T m H, 



y b m H± i 



y T m H± ; 



mi 



m 



1 + A^ CD) ) 



bt J j 



(2.6) 
(2.7) 
(2.8) 
(2.9) 

(2.10) 

(2.11) 



where N c = 3 is the color factor and we neglect the masses of the fermions other than the top 
quark. Here, we denote the QCD corrections to the individual decay processes as A^ CD ^ 
and A^ C . Importantly, the QCD corrections are calculable once the Higgs masses are 
determined |3*81 EH] ■ Thus, if we can measure the decay widths and branching ratios of the 
heavy Higgses, we can constrain the Yukawa coupling constants. In addition, the total decay 
width of the heavy Higgses are given by 



^A->bb + r^TT , 
^H^bb + ^H-^tt, 

^H±^bt + rg±_K 



(2.12) 
(2.13) 
(2.14) 



As we will see below, in the rapid-annihilation funnels, decay widths of heavy Higgses are 
relatively large. Thus, from the invariant- mass distribution of the decay products, the decay 
widths may be measured at the ILC. 

In summary, for the study of the heavy Higgs sector in the large tan/3 region, there 
are five relevant parameters so far: three masses rriA, ttih, fnn ± , and two Yukawa coupling 
constants y b and y T . As we will see, the relic density of the LSP can be also well calculated 
in the rapid-annihilation funnels once these parameters are fixed. In the following, we will 
discuss how well these parameters are experimentally determined. 

Although our method can be applied to various cases, the whole parameter space is 
too large to be studied. Thus, we choose several parameter points and see how well we can 
experimentally measure parameters at the ILC in those cases. For simplicity, we generate the 



may also decay into the gauge-boson or Higgs boson pair, but those processes are extremely suppressed 
in the large tan/3 region that we are interested in. 
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Point 1 Point 2 

~m Q 465.0 GeV 418.0 GeV 

mi/a 362.0 GeV 279.0 GeV 

A 60.0 GeV -12.0 GeV 

tan /3 50 48 

Hh 444.5 GeV 357.7 GeV 

m A 400.0 GeV 350.0 GeV 

m H 401.4 GeV 351.0 GeV 

m H ± 412.9 GeV 363.6 GeV 

T A 20.7 GeV 17.1 GeV 

T H ± 19.7 GeV 14.9 GeV 

B A _ bl 0.896 0.885 

m x o 144.4 GeV 109.6 GeV 

m x J 277.3 GeV 208.3 GeV 

m x o 450.4 GeV 364.8 GeV 

m x o 467.0 GeV 383.0 GeV 

m x ± 277.8 GeV 208.7 GeV 

m x ± 467.2 GeV 383.3 GeV 

tt L sph 2 0.113 0.113 

Table 1: Underlying parameters to be used in our study. 



mass spectrum of the MSSM particles by adopting the mSUGRA boundary condition which 
is parameterized by the following parameters: universal scalar mass mo, unified gaugino 
mass mi/2, tri-linear coupling coefficient A , tan/3, and the sign of the SUSY invariant Higgs 
mass sign(//#). Although we adopt the mSUGRA model as the underlying theory to fix the 
mass spectrum and mixing parameters, we consider a procedure to experimentally measure 
the parameters without assuming the mSUGRA. Since we are interested in the implication 
of such analysis to the study of the dark-matter physics, we adopt underlying parameters in 
the rapid-annihilation funnels. In Table ^ we list the underlying parameters which we use. 
(We also show the density parameter of the LSP fl^sp with those underlying parameters.) 



3 Heavy Higgses at the ILC 
3.1 Outline 

Now we discuss the study of the properties of the heavy Higgs bosons at the ILC. In partic- 
ular, we are primarily interested in how and how well we can constrain the five parameters 
listed in the previous section by combining the information from the neutral and charged 
Higgs production processes. For this purpose, we adopt several reasonable approximations 
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to simplify our analysis as we explain below, instead of studying the details of individual 
production and decay processes. 

As we mentioned in the previous section, decay rates of the heavy Higgses are sensitive 
to the Yukawa coupling constants. Thus our strategy is to study branching ratios and the 
energy distribution of the final-state fermions produced by the decay of the heavy Higgses. 
Since the (total) decay rates are much smaller than the masses of the heavy Higgses, we 
can also measure the masses of the heavy Higgses from the energy distributions of the decay 
products. 

Since we will try to determine the decay rates from the energy distribution of the final- 
state particles, we do not use the narrow width approximation for the heavy Higgs produc- 
tion cross section. In our analysis, the following type of the processes play a significant 
role: e + e~ — > <p\<p* 2 ~^ /1/1/2/2 where, in our case, (^1,^2) is or (H + ,H^) and /'s 

are quarks and leptons. Let us denote the invariant-mass squared of the faf! system as 
Sf-fi. Then, for the processes we are interested in, the differential cross sections are well 
approximated by the following formula* 4 

1/2 f, _ / 1/2 x 1/2 f, _ / 1/2 x 



s , 



2 



ds fiA ds h% 7x2 ( s hK - m li) 2 + T li m li ( s hf 2 ~ m ' 2 ) 2 + r ' 2 m ' 2 

+ <p 2 ), (3.1] 



where the second term should be omitted if /1 and f 2 are identical. Here, a(s; s^/*, Sf 2 p / 
is the cross section for the process e + e~ — > f\f* 2 with the masses of ip* and tp*, being set 
equal to sj/j, and s~j^, , respectively. In addition, f w -+/j/?( s /j?) denotes the decay rate of 
if* whose mass is identified as sYj>; f ^^.^(sV?,) can be obtained from Eqs. (j2.fi | — (12.11)) 

% l l/2 

by replacing the physical mass of the Higgs by s/ p . Explicit formulae of a will be given in 

JiJi 

the following subsections. Notice that, when T^. <C m ftl the total cross section becomes 



/ ds hf[ ds f2f! 2 ds _ ids - - ( B <pi->flfi B W-*hf 2 - + B ^flfi B ^f2f 2 ) a e+e-^ lV2 , (3-2) 
Ji /1 / 2 f 2 



where 



a e +e--*yxLp 2 {s) = ^(sjmj^mj,). (3.3) 



As suggested by Eq. 1)3.1)1 ■ the distributions of the invariant masses of the final-state 

1 /2 

fermions have peaks around sjp ~ m ft . In addition, the distributions become broader 



# 4 For the process we consider (i.e., e + e — > /1/1/2/2); there are in fact several diagrams where some of 
the fermions are not emitted from the virtual Higgs bosons (like e + e~ — > /*/(, followed by /j* — > f\+<f* and 
ip* — > fzfz). In the parameter region where sVl and s/^, are both close to the Higgs mass(es), however, 

Jl/l /2J2 

we have checked that the cross section is dominated by the process e + e~ — > — * /i/i/a/2- Thus, we use 
Eq. jOJ for simplicity. 
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as the decay rate T^. becomes larger. Thus, in the heavy Higgs production processes, the 
invariant-mass distributions of the final state fermions have important information about 
the masses and the decay widths of the heavy Higgses. 

It should be, however noted that the observed invariant-mass distributions are deformed 
by several effects. One reason is the energy loss by the neutrino emission. With the leptonic 
decays, some fraction of the initial energy of the bottom quark may be carried away by 
energetic neutrinos and the energy of the 6-jet becomes underestimated. Another important 
effect is from the resolution of the calorimeters (in particular, the hadronic calorimeter). In 
order to simulate these effects, we perform a Monte-Carlo (MC) analysis. Here, we first 
generate parton-level events for a given value of the collider energy y/s, which we take to 
be 1 TeV throughout this paper unless otherwise mentioned. In each event, in general, 
momenta of bottom quarks pf. , other lighter quarks and leptons are generated. If 
the final-state partons contain the top quark, its decay (as well as the subsequent decay of 
W /± ) is also taken into account at this state. Here and hereafter, the superscript (0) is for 
momenta of the primary partons before the hadronization and cascade decay. In order to 
consider the energy loss of the bottom quark in the leptonic decay events, we used ISAJET 
package [40 to follow the hadronization and decay chain of the bottom quark. For individual 
primary bottom quark with (initial) momentum p^\ we calculate the fraction of the visible 
energy after its hadronization and decay. Then, the 6-jet (after the hadronization and 
the decay) is treated as the jet with the momentum fbiPb^- In addition, because of the 
detector resolution, distribution of the observed energy of the jet is smeared. In order to 
determine the observed momentum of the individual jets, we adopt the energy resolution 
of the hadronic calorimeter to be o-^jE = 40 %/y/E (with E being the energy of the jet) 
(HUT]. As a result, with the generated momenta of the quarks pjj® and p^\ we determine 
the observed energy p^ and p qi . (We also denote corresponding observed energy as E^ and 
E qi .) 



3.2 Neutral Higgs production 

Let us first consider the neutral Higgs production e + e _ — > AH. In the large tan/? region, A 
and H both dominantly decay into bb and rf. Thus, there are three types of final states:* 5 



e+e" -> A*H* -> { bbrr . (3.4) 




The cross section for the process e + e L R —>■ A*H* is given by 



# 5 For simplicity, sometimes we do not distinguish particle and anti-particle when there is no confusion. 
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where gz = \jg\ + g\ with g\ and gi being gauge coupling constants for the U(l)y and 
SU(2)l gauge groups, respectively, mz is the Z-boson mass, a is the mixing angle in the 
CP-even Higgs sector, 



v AH = - 



i 2 2 \2 a 2 2 

v s — m A , — m^J — 4m A »m H » 



and, for the left- and right-polarized electron beam, 

, 92 - 9i , 

JL — , JR — 



9z 



9± 

9z 



(3.6) 



(3.7) 



We take sin 2 (/9 — a) = 1 since we are interested in the limit where the heavy Higgses are 
much heavier than the lightest Higgs. 

In our analysis, we use the first two types of events: 6666 (46) final state and bbrr 
final state. We do not consider tttt final state because the branching ratios of the heavy 
neutral Higgses into the rr final state are fairly small (~ 0.1) so the number of 4r events is 
suppressed. 

The branching ratios B A ^ bb and B H ^ bb are ~ 0.9 so sizable number of 46 event is ex- 
pected. At the ILC, such an event will be identified by the following features: 

• 4 6-tagged jets. 

• Small missing energy. 

• No isolated leptons. 

Once observed, the distribution of the invariant mass of two 6-jets provides important 
information about the properties of the heavy Higgs bosons. Energy and invariant-mass 
distributions of the primary partons are relatively easily calculated. As we mentioned in 
the previous subsection, however, observed distributions are different from the primary ones 
because of the leptonic decay of the bottom quark and also because of the resolution of the 
calorimeters. We are particularly interested in the distribution of the invariant mass of two 
6-jets which originate from (virtual) neutral Higgs bosons. In order to pair the 6-jets, we 
use the fact that the masses of the two (heavy) neutral Higgses are expected to be quite 
degenerate. Consequently, the total energy of 66 system from the virtual A becomes very 
close to that from the virtual H. Thus, for the event e + e~ — > bib'^b^ (where the 6-jets are 
ordered so that the observed energy be E bl < E b2 < E b > 2 < Ey ) , we define the following two 
invariant masses: 



= (Pbi + Pb[) 2 , rn 2 bb ,2 = (Pb 2 + Pb' 2 f 



(3.8) 



In general, we can perform statistical analysis based on the distributions of these two invari- 
ant masses. For simplicity, however, we use "averaged" distribution 



it, 



dm 



bb 



If 



dm, 



bb,i 



dN< 



46 



dm bb Am bb 2 



+ 



J rn bbi 2=m bb 



dm. 



66,2 



46 



_dm bbjl dm bb! 2 _ 



(3.9) 



m bb} i=m bb 
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where N 4b is the number of 4b events. As we will see, even from this averaged distribution, 
we obtain significant constraints on the masses and decay rates of the neutral Higgses. 

In order to determine the distribution given in Eq. ()3.9)1 . we calculate the "transfer 
functions" for the 4b events T4 b (m bb ; s blbl , s b2b2 ); T± b is the distribution of the observable rribb 
for the event with s blbl and s b2b2 being fixed. In calculating the transfer function, several 
kinematical constraints are taken into account to eliminate standard-model backgrounds. In 
the 4b events, all the beam energy is carried away by the bottom quarks. In addition, for 
the process e + e _ — > A*H* followed by A* — > bj)i and H* — > & 2 &2, E bl + E bi becomes close to 
E b2 + Ei 2 in most of the cases. The kinematical cuts based on these can be used to eliminate 
some of the backgrounds; here, for the event e + e" — > bib'^b^ with E bl < E b , 2 < E b > 2 < Ey , 
we adopt 

• 960 GeV < E bl + E b[ + E b2 + E b > 2 < 1040 GeV. 

• 470 GeV < E bl + E b[ < 530 GeV and 470 GeV < E b2 + Ey 2 < 530 GeV. 

• There is no leptonic activity (with energy greater than 25 GeV), in particular, in the 
6-jets. 

• Invariant masses of any two of the jets are larger than 130 GeV. 

With the excellent 6-tagging capability of the ILC, we expect that significant fraction 
of the 6-jets can be identified. Tagging of all four bottom quarks will help eliminating 
standard-model backgrounds. Using the transfer function, we obtain 

dN^b 4 r ,da e+e -^ A * H ,_> bbbb . . 

= e b L J ds blbl ds b2b2 T 4b {m bb ; s 6l5l , s b2 b 2 ) ^ _ ^ _ , (3.10) 

with L being the luminosity. Here, €b is the tagging efficiency of single bottom quark; we 
approximate that the 6-tagging efficiency is independent of the number of 6-jets and take 
e b = 0.7 in our numerical calculations PUI21IH]- Then, in the statistical analysis, we calculate 
the number of 4b events integrated over some intervals of the invariant mass: 

N d) = r- +1) dm J N *» 

4b Jm$ drribb 

where m bb and m^ 1 ^ are lower and upper bounds of the i-th bin. 

Although we imposed several kinematical cuts, there may still remain backgrounds. Since 
we required that the total energy of the final-state jets is close to y/s as well as that 4 b- 
tagged jets exist, we suppose that the dominant background is from the events of the type 
e + e~ — > bbbb. We use COMPHEP package [12] to generate such events and, using the 
kinematical cuts we discussed before, we estimate the number of backgrounds. We found 
that, however, the number of background from the process e + e~ — > bbbb is negligibly small 
(~ 0.01 event/5 GeV/1 ab _1 ). There may also exist another class of backgrounds which 
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(3.11) 



arise from the mis-identification of the 6-jet. Study of such background requires detailed 
study of the detector effects and we leave such a study for a future work.* 6 

In Fig. we show the number of events and backgrounds in each bin. Since A and H are 
quite degenerate, it is convenient to define the "averaged mass" and the "mass difference" 
as follows: 

rfiAH = ^(m A + m H ), Am AH = ^(m A - m H ). (3.12) 

(We will see that error in fa ah is much smaller than that in Am A H-) As one can see, the 
distribution of m bb is peaked at around m bb ~ rh AH . In addition, width and height of the 
distribution depend on the mass and decay parameters. Thus, from the observed distribution 
of m bb , mass and decay width of the neutral Higgses will be determined. 

Some fraction of the neutral Higgses also decay into rf pair, so we can also use the events 
with bbrr final state. Here we only use the hadronic decay mode of r to identify the r jets; 
when the r lepton decays hadronically, we obtain a jet with low multiplicity. Using this 
fact, we assume that the hadrons from r can be well identified and distinguished from the 
jets from the direct production of quarks. Then, the bbrr event is specified by the following 
features: 

• 2 6-tagged jets (b and b') with the total energy of the system being close to \\fs. 

• 2 energetic isolated jets with low multiplicity. 

Since, in the signal event, two 6-jets originate from single (virtual) neutral Higgs, we use 
the distribution of the invariant mass of two 6-jets to determine the masses and widths of 
the neutral Higgses: 

m lb = (Pb + Pb') 2 - (3.13) 

For the bbrr event, we also calculate the transfer function T bbrr with MC simulation. Then, 
we obtain 

dN bbrT 2 2 f . ,da e + e -^ A * H *^ bbTf 



dm 



e 2 b LB^ had f ds bb ds Tf T bbTT (m bb ; s blbl , s b2h ) ° e+e ^H^bbrr ^ ^ ^ 

bb •> as bb as T f 

where -B r ^had — 0.65 [21] is the hadronic branching ratio of r. In calculating T bbrT , we take 
account of the following kinematical cuts: 

• 470 GeV <E b + E v < 530 GeV. 

• In the 6-jets, there is no leptonic activity with energy higher than 25 GeV. 

• Angle between two jets from r leptons is larger than |7r. 



# 6 Invariant-mass distribution of the 46 event has been also studied in 221- where it is also shown that the 
background for the 46 event is well below the signal. 
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Figure 2: Number of 46 events in each bin (with the width of 5 GeV) with L = 1 ab~ . 
The histogram shows the result for the Point 1: rfiAH = 400.7 GeV, Aitiah = —0.7 GeV, 
Ta = 20.7 GeV, and B A _^ = 0.896. The short horizontal lines with "x" are those with 
one of the four parameters being changed: (a) rfiAH = 395 GeV, (b) Aitiah = 10 GeV, (c) 
T H = 15 GeV, and (d) B A _ bi = 0.85. 



For the bbrr event, the tt production process provides irreducible background; e + e~ — > tt 
followed by the r-leptonic decay of both the lV ± -bosons produced by the top decays. For 
this process, we generate the events and estimate the number of background of this type. We 
found that the number of events from this type of background is (0.1 — 1) event/5 GeV/1 ab _1 , 
which is below the number of signal event. As well as from the tt production process, we 
may also have backgrounds from the four-body production process e + e~ — > bbrf, which is 
also possible in the standard model. To see how large it is, we generate such event with 
COMPHEP package and estimated the number of background. We found that the number 
of such background events is very small after imposing the kinematical cuts, ~ 0.01 events 
per each bin (with L = 1 ab _1 ) which is well below the number of signal event. 
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Figure 3: Same as Fig. but for the bbrr events. The background from the tt production 
process is also shown in the darkly shaded histogram. 

Invariant-mass distribution of the bbrr event is shown in Fig. El As in the case of 4b 
event, the distribution is peaked at around the heavy neutral Higgs masses. The number of 
bbrr event is, however, proportional to B A _, b i(l — B A _^ b i) and hence is suppressed compared 
to that of the 4b event. As a result, 46 event is statistically more important than bbrr event. 

Using the above results, we can estimate the expected errors in the measured values 
of the physical quantities. For this purpose, we first determine the underlying values of 
the fundamental parameters. For the neutral Higgs production processes, there are four 
free parameters; here we use m A , m H (or m A H, ^ m AH), and B H ^ h i as parameters 
to specify the point in the parameter space. We use the parameter sets given in Table 
[U as underlying parameters unless otherwise stated. With the underlying parameters, we 
calculate the expected number of 4b and bbrr events in each bin, which we denote and 
Nbb TT - Then, in order to see how well the underlying values can be determined, we calculate 
the number of events in each bin, and N^ TTl for postulated values of m A , tuh, and 
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Figure 4: Contours of constant ^Xn on ^ m AH vs. tjiah planes; #Xn = 1, 2, and 4 (with 
L = 1 ab _1 ) from inside and the "star" on the figure indicates the underlying point. As the 



underlying point, we take the Point 1 in Tabled For the calculation of N± b and N^ v in 
Eq. ([3. 15)1 . we take Ta and B A _> b i to be (a) Ta = 20.7 GeV and B A ^ bb = 0.896 (underlying 
values), (b) T A = 19.7 GeV and B A ^ bb = 0.896, (c) T A = 20.7 GeV and B A ^ bb = 0.88, (d) 
T A = 19.7 GeV and B A ^ bb = 0.88. 



Bu^bli an d see if N^ b and N bbTT are statistically consistent with Njjjjj and N bbTT . Here, we 
define the 5x 2 variable for the neutral Higgs production as 



E 



W 

16 



N, 



«N2 



16 





bbrr 



A7-(«) \2 
iV 66rrJ 



46,BG 



4- N 



(3.15) 



66rr.BG 



where N^ bBG and N bbrT ^ BG are numbers of backgrounds for the 46 and bbrr events, respec- 
tively. 

With 5x^5 we estimate how well we can determine the physical quantities. Since 
depends on four independent parameters, we fix two of them and show the behaviors of <5xn 
on several two-dimensional hyperplanes. First, we consider the constraints on the masses of 
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Figure 5: Contours of constant on AmAH vs. Ta planes; = 1, 2, and 4 (with 
L — 1 atT 1 ) from inside. We fix ttiah and B A ^bl to be (a) m^n = 400.7 GeV and B A ^bl = 
0.896 (the underlying values), and (b) ttiah = 400.7 GeV and B A _+ h i = 0.88. For the 
upper contours, the underlying values are given by the Point 1. For the lower contours, the 
underlying value of Ta is taken to be 10 GeV. (In (b), there is no countour for 5xn = 1 for 
the case of Ta = 10 GeV.) 

the neutral Higgses. For this purpose, we consider the hyperplane with Ta and B A -^bl being 
fixed and show the contours of constant <?>Xn on Am ah vs. ttiah plane. In Fig. HI we show 
the results. In particular, in Fig. Ufa), we show the behavior of <5xn ° n the hyperplane with 

and B A -^bl being equal to their underlying values. As one can see, if we consider the 
contour of <5xn = 1> error of rfiAH is ~ 1 GeV while that of AmAH is ~ 6 GeV. We also 
show the contours of constant on other plane. We can see that the contours for = 2 
and 4 are almost unchanged while the region with Sxh < 1 is reduced. From the fact that 
the error in ttiah is significantly smaller than that in Aitiah-, we expect that the error in the 
determination of the masses of heavy neutral Higgses is dominantly from that of Am A H- 

Next, let us consider the uncertainties in other quantities. In Fig.0 we show the contours 
of constant <$Xn on the hyperplane of friAH an d B A ^bl being fixed. From this figure, we can 
see that the uncertainty in T a is 2 — 3 GeV when the underlying value is T a — 20 GeV. In 
fact, the uncertainties are sensitive to the underlying value of T a- To see this, we also show 
the results for the case where the underlying value of T a is 10 GeV (with other underlying 
parameters being unchanged). In this case, we can see that the errors in Am ah and T a are 
reduced. 

Using the variable, we can also constrain the branching ratios. In Fig. we plot the 
contours of constant <5x^ on Amn vs. B A ^bb plane. We can see that B A _> b i can be constrained 
at the level of ~ 5 % and hence we can conclude that B A ^bl c& n be well determined. 
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Figure 6: Contours of constant for the Point 1 on AniAH vs. B A _+ h i planes; <5xn = 1> 2, 
and 4 (with L = 1 ab™ 1 ) from inside. Here, ffiAH and Ta are fixed to be (a) in ah = 0.7 GeV 
and Ta = 20.7 GeV (the underlying values), and (b) ffiAH = 0.7 GeV and Ta = 19.7 GeV. 



3.3 Charged Higgs production 

Next, we consider the charged Higgs events. The useful decay modes in the large tan/3 case 
are H — > bt and tv t . Thus, in the charged Higgs events, relevant final states are 

( btbt 

e + e~ -> H + *H-* -> I btrv T . (3.16) 

[ tv t tv t 

We expect significant amount of irreducible backgrounds for the tv t tv t event in particular 
from the W + W~ production. Thus, we do not consider this mode. 

First, we consider the btrv T final state. After the decay of the top quark, we obtain W ± - 
boson and b. In order to determine the mass of the charged Higgs, we only use the hadronic 
decay of the W^-boson. In addition, for r, we again use only the hadronic decay mode to 
identify the r-lepton events. Then, the relevant btrv T event has the following features: 

• 2 6-tagged jets {b\ and 62)- 

• 2 non-Mike jets (qi and q 2 ). 

• 1 energetic isolated jet with low multiplicity. 

To eliminate backgrounds, we also impose several kinematical constraints. Since the 
event is dominated by the back-to-back production of the charged Higgs bosons (which are 
almost on-shell), total energy of the bt system is close to hy/s in most of the signal events. 
In addition, if we consider the system consisting of q±, g 2 and one of the 6-jet, its invariant 
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mass becomes close to m t . (We may be able to use the fact that the invariant mass of the 
system consisting of qi and q 2 system is close to m w although, in our numerical analysis, we 
do not take into account such constraint.) 

In the btrv T event, all the hadronic activities except the isolated hadrons from r are from 
one of the charged Higgs, so the invariant mass of such system contain information about 
the mass and decay width of H^. Thus, we define 

m bt = (Pfci + Pb2 + Pqi + Pq 2 ) 2 > ( 3 - 17 ) 

and calculate the distribution of mu by MC analysis. In deriving the distribution of this 
invariant mass, we impose the following cuts: 

• 470 GeV < E bl + E b2 + E qi + E q2 < 530 GeV. 

• 150 GeV < m biqiq2 < 200 GeV or 150 GeV < m b2qiq2 < 200 GeV, where m 2 bmq2 = 
(Ph +P Ql +Pg 2 ) 2 - 

• No leptonic activity (with energy larger than 25 GeV) in the 6-jets. 
Then, we calculate the transfer function T btTUr to estimate 



dN btTUr 2 r •,da e + e -^ H +*H~*-*btru T /„ 10 n 

— i — e b L1 ^T^ha,d^W ± ^qq I dS b t<lS TUT 1 btTU T \ m bt'i $bt, S TV ) , ^d.iSJ 

dm bt J dsuds^ 



where Bw±^ q q — 0.676 [21] is the branching ratio of the hadronic decay of W -boson. Here 
the cross section for e + e~ — > H + *H~* is given by 

3 / 2 „2 _2 - \ 2 



sv 6 H ± ( e z gj - gj 1 

TS — r J L > R 

487T V s g z s-m z 



° e +e7 „-*H+*H-*( S 'i m H+*, m H-*) = I T + ^ " ~ ' ( 3 " 19 ) 



where e is the electric charge, and Vh± is obtained from Eq. (|3.6|) by replacing uia* — > m H +* 
and triH* — * van-*. 

In Fig. [7[ we plot the distribution of the invariant mass m bt for several cases. First 
of all, as one can see, the distribution function is peaked at around ~ m.#±. In addition, 
the distribution becomes broader as the decay width becomes larger. Notice that the 
number of btTv T event is approximately proportional to B H ±^ bt {l — B H ±^ bt ). 

For the btrv T process, there exists irreducible background from the ti pair production 
process, if one of the boson decays hadronically while the other decays into r and v T . 
We have estimated the number of the background from ti pair production, generating such 
events. The result is also shown in Fig. [7| the number of background is well below the 
number of signal event. We also calculated the background from the standard-model process 
of e + e~ — > btTu T with the COMPHEP package with eliminating the contribution from the 
process e + e~ — > ti, we found that the number of background event of this type is extremely 
small. 
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Figure 7: Number of btrv T events in each bin. The shaded histogram shows the result for the 
Point 1: m H ± = 412.9 GeV, T H ± = 19.7 GeV, and B H ±^ bt = 0.885. The short horizontal 
lines with "x" are those with one of the parameters being changed: (a) m#± = 405 GeV, 
and (b) T H ± = 15 GeV. The background from the tt production process is also shown in 
the darkly shaded histogram. 



As a charged Higgs production event, btbt final state event is also available. In this case, 
we expect four 6-tagged jets and other hadronic and/or leptonic activities in the final state 
(from the decay of the W ± bosons). In the btbt events, however, the reconstruction of the 
invariant mass of the bt system of one side seems challenging because of the combination 
error or overlap of the jets (although it may be possible with some careful analysis). Thus, 
in this paper, we only use the total number of the btbt event to constrain the branching ratio, 
using the fact that the number of the btbt event is approximately proportional to B 2 H± ^ bt . 
Here, we estimate the number of btbt event as 

Nbtbt = 4>B 2 H±^btLc> e+ e- H- ■ (3.20) 

As in the case of the neutral Higgs, we calculate the Sx 2 variable; for the btrv T event, we 
calculate the number of events falling into each bin, which is classified by mu- The number 
of event in the z-th bin N b ^ Tll is calculated as a function of independent parameters (in the 
charged Higgs case, m H ±, T H ±, and B H ±^ bt ). For the btbt event, we use only the total 
number of events N btbt . Then, we define 

x 2 _sr ^btrur ~ N btTu T ) 2 , (N b tbt - N bt bt) 2 /„ 01 \ 



NiL + Nbf^a Nbtbt + N 



btbt, BG 



where the second terms in the denominators are the numbers of background events. 

In Fig. |Hl we show the contours of constant 5xc on m H ± y s- planes. We can see 
that, from the contour of 5xc — 1) the uncertainty in mn± is 1 — 2 GeV while that in T H ± is 



18 



30 
25 



% 20 



x 



15 



10 - 



30 



25 



I 20 

-H 

U 15 



10 




410 



415 

m H ± (GeV) 



410 



415 

m H ± (GeV) 



420 



Figure 8: Contours of constant 5xc on m H± vs. planes; Sxc 
L = 1 ab" 1 ) from inside. The branching ratio is fixed to be Bjj±. 



1, 2, and 4 (with 
= 0.885 for (a) and 

(c), and B H ±^ bt = 0.87 for (b) and (d). Upper two figures ((a) and (b)) are with tbtb mode 
while the lower ones ((c) and (d)) are without tbtb mode. 

~ 3 GeV. We can also see that the inclusion of the tbtb mode can help reducing the allowed 
region on the vs. r#± plane. In Fig. El we also show the contours of constant 5xc on 
rriH± vs. Bh±->h planes. We can see that the branching ratio Bh±^m is also well constrained 
by the study of the charged Higgs production processes. 

3.4 Combined results 

Now, we combine the results obtained in the previous subsections. For this purpose, we 
define the total 5\ 2 variable as 



^Xtot 



(3.22) 
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Figure 9: Contours of constant Sxc on m H± vs. B H ±^ bt planes; 5\% = 1, 2, and 4 (with 
L — 1 atT 1 ) from inside. Here, we take T H ± = 19.7 GeV for (a) and (c), and T H ± = 17.0 GeV 
for (b) and (d). Upper two figures are with tbtb mode while the lower ones are without tbtb 
mode. 

This quantity depends on five parameters (as well as on the underlying parameters); we use 
ttia, m H, "mH±, Ta, and B A _^ h i as free parameters. Once these parameters are fixed, all the 
Yukawa coupling constants (as well as the QCD correction factors) are determined in the 
large tan/3 region and the number of the neutral and charged Higgs production events can 
be calculated. As we saw in the previous subsections, averaged mass of the neutral Higgses 
ffiAH, the charged Higgs mass m#± and the branching ratios are relatively well determined. 
Thus, in this subsection, we concentrate on the uncertainties in the remaining two parameters 
which are most important for the calculation of the relic density of the LSP, i.e., Am ah and 

In Fig. Ell we plot the contours of constant Sxtot f° r Point 1 on Atuah v s- plane, 
on which rfiAH, m u ± i an d B A ^ h i are fixed to be the underlying values. As one can see, 
constraint on the ArriAH vs. Ta plane becomes more stringent as we add the information 
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(with L = 1 atT 1 ) from inside, ffiAH, t^h±i an d B A ^ b i are fixed to be the underlying values. 
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Figure 11: Same as Fig. but for the Point 2. 
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from the charged Higgs production events. Compared with Fig. we can see that the 
errors in Am^ and Ta are reduced by including the information from the charged Higgs 
productions. We also performed the same study for the Point 2, and the result is given in 
Fig. El 

Before closing this section, let us comment on the determination of the tan f3 parameter. 
Since the decay widths of the heavy Higgses are sensitive to tan (3, we can obtain information 
about tan (3 from the decay widths of the heavy Higgses P . In order to derive the value of 
tan f3, however, we need to relate the Yukawa coupling constants to the observed masses of 
the (third-generation) fermions. For the bottom quark, this is difficult since the radiative 
correction from the supersymmetric diagrams may significantly affect the bottom quark 
mass. Expecting that the radiative correction to the mass of r is small enough, however, we 
can use the relation 



tan/3 = — — . 



(3.23) 



Then, using the decay rate of the CP-odd Higgs boson, for example, we obtain 



tan/? 



'16^(1-5 



A-+bb) 



m A 



(3.24) 



and hence 



S tan (3 
tan (3 



A 



+ 



'5mA 



+ 



5B 



A^bb 



B 



A^bb, 



1/2 



(3.25) 



As we have seen in the previous subsections, uncertainties in Ta and Ba^tt = 1 — B A ^bl 
are both 10 — 20 % level, while that in uia is smaller (~ 1 %). Thus, in this case, the tan/3 
parameter can be determined with the accuracy of 10 % or so. 



4 Reconstruction of Dark Matter Density 

Now, we are at the position to consider the possibility of reconstructing the dark matter 
density in the rapid-annihilation funnels. In the parameter region we are interested in, pair 
annihilation of the lightest neutralino (which is the LSP) is dominated by the process with 
the s-channel exchange of the pseudo-scalar Higgs boson A. The cross section for this process 
is calculated as* 7 



^Va^-v^A s 



a " rn , a „ n — m 2 .\^ 4- Y 2 .m 2 . v ' 



m A [s x o x o - m A r + 1 A m : 



# 7 Strictly speaking, in Eq. (|4.1|l . we have to take account of the QCD correction to the decay rate when 

1/2 

s o o differs from tua- Such a QCD correction is calculable and, in this paper, we neglect such a correction. 
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where s x o x o is the invariant mass of two LSPs in the initial state, and UaxIxI * s ^ ne coupling 
between A and the lightest neutralino, which is given by 

y Ax o x o = l - (g x [U x o] u - g 2 [U x o] 12 ) (sin p [U x o] 13 - cos (3 [U x o] u ) . (4.2) 

Here, LLo is the unitary matrix which diagonalizes the neutralino mass matrix which is given 
by 



Mo 



( -m G i -givcosfi givsm.fi ^ 

— mc2 cos /3 — g>2^ sin (3 

—givcosj3 g2Vcos(3 \in 

\ g\v sin /3 — g2V sin /3 /x# / 



(4.3) 



with mci, rn.G2, and yiiij being the gaugino masses for £/(l)y and SU(2)l gauge groups and 
SUSY invariant Higgs mass, respectively. Thus, in the rapid-annihilation funnels, dominant 
contribution to the pair annihilation cross section of the lightest neutralino is calculated once 
rriA and T A as well as Uax°x are determined. In the previous section, we have seen that, by 
the study of the heavy Higgses, tua and T A are constrained fairly well. 

For the determination of VaxIxV we ne ed the mixing matrix U x o (as well as tan/5) which 
is calculated from the neutralino mass matrix. Importantly, the neutralino mass matrix de- 
pends on the gaugino masses and the SUSY invariant Higgs mass which can be determined 
once the chargino and neutralino masses are experimentally measured. At the ILC, measure- 
ments of these masses can be performed if the charginos and neutralinos are kinematically 
accessible. In particular, by studying the production processes at the threshold region (i.e., 
by the threshold scan), some of the masses of the charginos and neutralinos can be deter- 
mined with the accuracy of ~ 50 MeV jSj- In addition, from the kinematics of the decay 
products of the chargino and neutralinos (as well as sfermions, if kinematically accessible), 
mass of the LSP (i.e., the lightest neutralino) is also determined with 5m x o ~ 50 MeV [3]. 
Thus, if all the charginos and neutralinos are seen at the ILC, it will give us enough infor- 
mation to precisely determine Uax^xI- ^ some of the charginos or the neutralinos are too 
heavy to be experimentally produced, it becomes rather difficult to determine the neutralino 
mixing parameters. In such a case, we may have to perform some global fit using all the 
masses of the superparticles. Of course, some information from the LHC may be also used. 
Since our main concern is to study the properties of the heavy Higgses, we do not consider 
the detail of the parameters in the neutralino sector. Instead, we assume that, at the ILC, 
all the charginos and neutralinos can be produced and that their masses can be precisely 
measured. If the masses of all the charginos and neutralinos are measured with the accuracy 
of ~ 50 MeV, the uncertainty of VaxIx" ^ s ex P ec ted to be 0(0.01 %) so that the dominant 
uncertainty in the reconstructed value of ^lsp is from uia and Ya- 

To see how well we can reconstruct ^lsp, in Figs. H21 and [T31 we plot theoretically 
calculated value of ^lsp on rriA vs. Ta plane. (In our numerical calculation, we use the 
DarkSUSY package [121 to calculate f2 L sp-) Here, we fixed Uax^xI anc ^ m x? as their underlying 
values and calculated ^lsp as a function of tra and Ya- On the same figure, we also show the 
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for the Point 1. The numbers in the figure are the value of Q^sph 2 . Contours surrounding 
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constraint on the itla vs. Ya plane expected from the ILC. Here, assuming that the dominant 
error of is from Am ah, we neglected the uncertainty of ttiah- From these figures, we 
can see that the uncertainty of ^lsp is 10 — 20 %, which is comparable to the uncertainty 
of the dark-matter density determined from the WMAP data. 

We have not discussed possible errors of ^lsp originating from parameters other than m A , 
Ta, and y Ax o x o. In order to calculate the total pair annihilation cross section of the LSP, we 
also have to take account of other processes like the s-channel H° exchange and t-channel 
sfermion exchange processes. In the rapid-annihilation funnels, however, effects of these 
processes are subdominant. (We have checked that processes other than the CP-odd Higgs 
exchange change the value of ^lsp at most a few % in Points 1 and 2.) In addition, more 
importantly, once the properties of the heavy Higgses and the sfermions are determined at 
the ILC, cross sections for these subdominant processes are also calculated. Thus, we expect 
that the uncertainty in ^lsp from the processes other than the CP-odd Higgs exchange is 
less than a few % and that the dominant errors are from ijia and T^- 

5 Conclusions and Discussion 

We have considered the capability of the e + e~ linear collider for studying the properties of 
the heavy Higgs bosons in the supersymmetric standard model at the ILC. We concentrated 
on the large tan/5 region which is motivated, in particular, by explaining the dark-matter 
density of the universe (i.e., so-called "rapid-annihilation funnels"). We perform a systematic 
analysis to estimate expected uncertainties in the determination of the masses and widths 
of the heavy Higgs bosons. 

With the study of the invariant-mass distributions of the jets in the final state, we 
have seen that the masses, widths, and the branching ratios of the heavy Higgses are well 
constrained. Compared to the averaged value rfiAH = \{ m A + ttih), mass difference of the 
neutral Higgses, Am ah = ^(ttia — ttih), is more difficult to measure. Consequently, if we try 
to experimentally determine the masses of the neutral heavy Higgses, uncertainty in Am ah 
becomes the significant source of the error. 

In this paper, our primary purpose was to point out the strategy for the systematic 
study of the heavy Higgs bosons at the ILC and to estimate the expected uncertainties 
in the measurements of their masses and widths. Thus, we assumed that the 6-tagging 
efficiency and the energy resolution of the detector is well understood and we did not consider 
systematic errors from these. In addition, for the background, we take account only of the 
dominant physics background. When the study suggested in this paper will be performed, 
these points should be studied in more detail; in order to realize the detailed and precise 
study of the heavy Higgs bosons, the following will be necessary: 

• Good understanding and high efficiency of the 6-tagging. 

• High resolution of, in particular, hadron calorimeter. 

• Good understanding of the backgrounds. 
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We have also discussed the implication of the study of the heavy Higgses to the calcu- 
lation of the relic density of the LSP. In the rapid-annihilation funnels, the dominant pair 
annihilation process of the LSP is through the s-channel exchange of the CP-odd Higgs 
boson. In this case, the mass and width of the CP-odd Higgs boson should be determined 
for the precise calculation of the relic density of the LSP. We have seen that, if the LSP 
dark matter is realized in the rapid-annihilation funnels, the the dark matter density can 
be reconstructed with a very good accuracy of 10 — 20 % (see Figs. H21 and IT5|) . If the 
uncertainty in the mass or the width of the CP-odd Higgs can be reduced, we will have a 
better determination of I^lsp- In particular, the dominant source of the uncertainty in itla is 
from the determination of the mass difference of two neutral heavy Higgses. If the radiative 
correction to the Higgs potential is well understood, then the mass difference Attlah may 
be theoretically calculated. In this case, uncertainty in becomes smaller and a better 
determination of ^lsp is expected. In addition, we have also assumed that the radiative cor- 
rections to the pair annihilation processes of the LSP will become well studied by the time 
when the superparticles as well as the heavy Higgses are produced at the ILC. Although 
most of the radiative corrections are calculable, many of them have not been calculated yet. 
In order for the precise theoretical calculation of ^lsp? such a study will be very important. 

In conclusion, we have seen that the ILC is very useful not only for studying the properties 
of the new particles in particle-physics models beyond the standard model (in this case, 
the heavy Higgse in the sypersymmetric models) but also to have deeper insights into the 
evolution of the universe. In particular, the ILC may help answering one of the most serious 
misteries in cosmology, the origin of the dark matter of the universe. If the reconstruction of 
the dark matter density will be successful, it will provide us a better understanding of our 
universe up to the temperature of O(10 GeV). 
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